Theory of integer quantum Hall effect in bilayer graphene near neutrality point 
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A Hartree-Fock variational ground state for dual-gated bilayer graphene (BLG) in the presence of 
magnetic fields is proposed. We argue that near the charge neutrality point (CNP), a Neel as well as 
uniform magnetization can coexist in the presence of the magnetic fields. An infinitesimal Zeeman 
coupling restricts the Neel order to the easy plane and develops magnetization in the direction 
of the field. Many-body effects further enhance the size of the magnetization. Such interaction 
induced gap near the CNP displays a quadratic dependence with the field, when it is weak. It then 
crosses over to a linear one at intermediate strength of the field, and possibly persisting even at 
stronger fields. These features are in qualitative accordance with a recent experimental observation. 
Furthermore, with realistic strengths of the interactions, we establish an quantitative agreement with 
experiment for T < _B < 4.5 T. Upon tilting the field, the Neel order gradually vanishes and a pure 
ferromagnetic order emerges beyond a critical strength of the total field, which otherwise depends 
monotonically on its perpendicular component and the zero field gap. Formation of additional 
incompressible Hall states and the nature of the broken symmetry phases at other fillings, e.g., 
v — ±2, are discussed as well. 
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Introduction: Nature of the underlying broken sym- 
metries in quantum Hall systems has recently been an 
active field of research following the successful realiza- 
tion of the carbon-based layered materials, for example 
mono-layer graphene (MLG) and BLG. They respectively 
discern quantized plateaus of Hall conductivity at fillings 
1/ = ±(4n -1-1) and ±(4n -|- 2), at relatively weak mag- 
netic fields [y . The four fold degeneracy of each Landau 
levels (LLs) arises from the spin and the valley degrees of 
freedom. The additional orbital degeneracy of the zeroth 
LL in BLG, bears the signature of the quadratic quasi- 
particle dispersion around the Dirac points, opposed to 
the hnear ones in MLG @, @1- The distinct topologies 
of such simple band structures leave profound impacts 
on the role of electron-electron interactions in those sys- 
tems. Weak electron-electron interactions are irrelevant 
in MLG [j|, while they can cause various instabilities in 
BLG . The nature of the broken symmetric phases in 
BLG, is however sensitive to the range of interactions. 
A weak coupling renormalization group calculation pre- 
dicts that BLG can be found either in metallic nematic 
or insulating antiferromagnet (AF) phases, when inter- 
action is respectively long and short-ranged [y, lZ[. In 
recent experiments, the nematic |8i] as well as an insulat- 
ing phase has been observed, in the presence of single 
and dual gate, respectively. For simple onsite Hubbard 
model, the most dominant instability is towards the for- 
mation of a st agg ered pattern of spin among the layers, 
the AF phase p^ . 

In a half-filled bipartite lattice, for example BLG, the 
Stoner instability towards the formation of a ferromag- 
net (FM) order is preempted by the appearance of an 
AF phase. Therefore the double-gated BLG is gapped, 
possibly with an AF order. However, in the presence of 
magnetic fields, even an infinitesimal Zeeman coupling 
can modify the nature of the ordered phase. It projects 



the Neel order onto an easy plane, and concomitantly 
favors a magnetization in direction of the field, which 
unless otherwise mentioned is set to be perpendicular to 
the BLG plane. Many body effect further enhances the 
size of the magnetization. A similar ground state has 
also been proposed recently for neutral MLG subject to 
quantizing magnetic fields [ll|. A recent experiment [y], 
discerns a gap (2Ao) ~ 1 meV in dual-gated BLG in the 
absence of the field at the CNP {v = 0), which otherwise 
increases monotonically with the magnetic field (B), and 
closely resembles to the form 



Eqap — Ao 



A^ 



a^B^ 



(1) 



where a = 5.5 meV/T. Therefore, at low fields Egap gap 
varies quadratically with the field, which then reverts to 
a linear one at stronger fields. Dependence of the gap 
in the presence of on weak electric field, applied perpen- 
dicular to the BLG, rules out the possibility of a layer 
polarized state, corresponding to an imbalance of aver- 
age electronic density among the layers [IJ, U^- Here 
we analyze the scaling behavior of the above proposed 
ordered phase with the field. We show that developing 
a sizable magnetization in the direction of the field, a 
quantitative agreement with the experiment can be es- 
tablished, at least when the field is weak. 

The ground state: Hereafter we assume a uniform 
background of electronic density and Neel (N) order as 
well as magnetization (M) in the BLG system. The 
ground state energy per unit area of the BLG at half- 
filling can therefore be written as [Ijl 



N' 



M' 



Egr — h h Eq 



N,M 



(2) 



-Bo N, M is the ground state energy per unit area of 



the effective single-particle Hamiltonian 

Hhf = /2 «) i^o - (^ • ^) (g) 70 + (A + M) CTa (^ h, (3) 
where 



Ho = 72 



2m* 



71 



2m* 



(4) 



m* sa 0.033me is the effective mass of the parabolic dis- 
persion around two inequivalent Dirac points iK, with 
K = (l, l/Vs) (27r/a-\/3) @- rne is the mass of ordinary 
electrons and a is the lattice spacing of the underlying 



honeycomb lattice in each layers, ttj = (- 



idj 



Aj ) , with 



j — x,y and we set the magnetic field (B) to be perpen- 
dicular to the BLG plane, i.e., B = esijdiAj. Therefore 
the magnetization is kept only in the direction of the field, 
M = Mz. A = Xz^c is the Zeeman coupling, where Wc is 
the cyclotron frequency and Xz ~ 0.015 for BLG. Hhf is 
written in a rotationally invariant basis ^ 
where ^I = 



(*T'*i) 



Vl,cr{K + cf), V2,a{K + q) , Vi^ai^K + q) , V2.,cr{-K + q) 



(T =t, I corresponds to the projections of spin along the 
z-direction. « = 1,2 corresponds to two layers of BLG. 
The mutually anti-commuting four component Hermitian 
gamma matrices are 70 — <jq ® a^, ll-,~ '^3 ® '^2,12 = 
o'o <8i CTi , 73 = CTi (g) 0-2 , 75 = 0-2 (8) CT2 [13 ■ The spectrum 
of Hhf is therefore comprised of a set of Landau levels 
at well separated -ilen.(J^ where 



Ni 



1/2 

n{n - 1)ujI + N'(^ + ct (A + M) 



1/2 



(5) 

with N±_ = [Ni,N2) ,N'i = N\\ and cr = ±1, with degen- 
eracies per unit area I/tt^^ for n = 2, 3, 4, • • • and l/27rZ^ 
for n = 0, 1 [lg|- Ib — yific/eB)^ is the magnetic length. 
At half-filling, all the states at negative (positive) ener- 
gies are filled (empty). Therefore, 



E. 



N,M 




2ttII 



eia 



2^e,, 



(6) 
To find the configuration of TV to minimize the Hartree- 



Fock ground state energy, we choose 



pendent variables. Then OEq 



N,M 



N\ and iV|| as inde- 
/dNii = yields 



Y.<y{X + M) 



a=± 



E 

n=0,l 



^-E 



2iVii 



i>2 e 



no^^+n{n-l)u 



(7) 



The left hand side of this equation is negative definite 
function of N\\ and vanishes only for A^|| = 0. Therefore, 
in the presence of even infinitesimal Zeeman coupling the 
Neel order is projected onto the easy plane (A^n — 0), per- 
pendicular to the direction of the magnetic field. This 
configuration also corresponds to the minima of the en- 
ergy. An identical ground state has also been proposed 
for MLG previously jllj. Therefore, in the presence of 
the field, the total gap at charge-neutrality point [v = 0) 
is 



Egap{y = 0) = ^Nl + (A + Mf 



(8) 



Next we find the self consistent solution of \N±] and M 
as a function of the magnetic fields. 

Gap equations: One can now arrive at the self consis- 
tent gap equations by minimizing Egr w.r.t the Neel and 
the magnetization orders, separately. After some tedious 
otherwise straightforward algebra, in terms of a new set 
of variables 5^ = A^_i/Ao,(5„i = M/Aq,/? = Wc/Aq, the 
minimization of Egr w.r.t. N±_ and m respectively yield, 
either 5n = ^ or 

2 '-'^ 

9A 



a=±l 



^di 



+ 



2/32 



yiX(e,^iV,J™,A,/3,CT) ' ^^+fi2y2 



(9) 



and (5„, = or 



2<5™/3-i 



9F 



E 

cr=±l 



[ae+vOd^ 



22//32 



„ (10) 

takmg 8gA,F/2Tr -> gA,F [Hl- Here y = Xz + S„i/3 \ and 



X{t6N,P,y) = 




(11) 



In last two equations we neglect the contributions, not 
relevant at low fields, i.e., /? ^ 1. As the ultra-violet 
cut-off 51 ^- 00 the right-hand side of Eq. [5] exhibits 
a logarithmic divergence, whereas the next equation is 
completely divergence free. Such an outcome is directly 
related to the LL spectrum in Eq. [5] . The AF order and 
magnetization, respectively shift and split all the LLs 
at finite energies, whereas both develop gap within the 
zeroth LL. A similar spectrum of the LLs is found in 
MLG [ll|. To ensure the cut-off independence of the 



order parameters, we substitute 



1 
9A 



d^ 



ve+^o 






de 



v^ 



1 



(12) 



in Eq. [3] , where Aq is the zero field gap of the Neel or- 
der. On the other hand, gp remains as a free parameter. 
To evaluate the right-hand-side of the gap equations we 
reformulate them by using the Feynman parametrization 

M 



1 



dy 



y/y{i-y) 



1 



yA+il-y)B^ 



(13) 



One can then perform the integrals by using the "Eu- 
ler substitution" ^19.] . In this prescription, we manage 
to evaluate all the divergent integrals analytically. The 
remaining finite integrals are computed numerically. 
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FIG. 1: Left: Magnetization {5m = M/Ao) in neutral BLG 
for QF = 0.58 (red), 0.6 (black), 0.63 (blue), 0.65 (green). 
Right: AF (5n = A^x/A,,) order for gp = 0.58. Here /3 = 

LJc/Ao 

Numerical strategy: To find the self-consistent solu- 
tions of two order parameters, we here pursue the follow- 
ing numerical approach. For each field strength, first we 
find the zeros of the Eq. [5] in the {5N,5m) plane, yielding 
a closed curve. Taking the ansatz 5n — tan0 5m, corre- 
sponding to a family of rays passing through the origin 
otherwise enjoying intersections with that closed curve, 
we compute only positive roots ( i.e., < < 7r/2 ) of 
Eq. [9] . Plugging these roots in Eq. [10] , we search for the 
self-consistent solutions of 5^ and 5„n for various values 
oi QF- 

Scaling: The AF order is finite at zero field and grows 
quadratically with the field, if the field is weak (/3 ^ 1). 
The curvature of the quadratic variation of the AF order 
decreases with increasing gp^ though not considerably. 
The magnetization, on the other hand, varies linearly 
with the field, but vanishes in its absence. As gp in- 
creases, magnetization is enhanced. This features are 
depicted in Fig. [1]. Hence to the leading order, the gap 
near the CNF Egap{v — 0) scales quadratically with the 
field, when it is weak, as expected from Eq. [5]. The scal- 
ing of the total gap in the weak field limit for gp — 0.6 
is shown in left column of Fig. [2]. A similar qualitative 
picture has been attained previously considering only an 
AF order. However, the curvature is four times smaller 
than the observed one [y|. As gp gets stronger, the cur- 
vature of the scaling function increases, shown in right 
column of Fig. [5] . The coefficient of the quadratic term 
are found to be 0.13, 0.17, 0.25, 0.45 respectively for 
gp = 0, 0.6, 0.63, 0.65 from numerical fitting. For smaller 
gp{< 0.5) the curvature does not change considerably 
from the one with a pure AF order, i.e., gp = 0. To 
establish the connection of our result with the recent ex- 
periment 9], one must take into account that the zero 
field gap in Eq. [T] is 2Ao [9], while in our analysis it 
reads as Aq. The coefficient in front of /3^ in the scaling 
of interaction induced gap at z^ = 0, is then found to be 
0.5 in experiment, while that for gp = 0.65 gives 0.45 
from a numerical fitting, see Fig. [5]. Therefore by de- 
veloping a finite magnetization, ^ (10"'^ — 10~^)^b per 



electron when i? '-^ 1 T, where fj,B is Bohr magneton, the 
quadratic scaling of the interaction induced gap in weak 
fields, is found in reasonable agreement with the experi- 
ment •9,2(f\. On the other hand, at intermediate strength 
of the field (/? > 0.6 for gp — 0.65) the gap starts to de- 
viate from the quadratic variation, and enters the regime 
of a linear variation with the field, shown in Fig. [5]. 
The field strength associated with such crossover reduces 
with increasing gp . However, the slope of the linear scal- 
ing function remains close to the one found previously for 
a pure AF state. To be precise, the slope of the scaling 
function increases slightly, though it is still off by a factor 
of 1.85 from the observed one. The same quantity for a 
pure AF gap is found to be off by a factor 2 Q . 

The role of the single-particle Zeeman coupling, even 
though small, is crucial. We also compute the self consis- 
tent solution of Sn and (5,„, setting A = 0. A non-trivial 
magnetization then develops only for gp > 3. The Stoner 
instability from N± 7^ 0, rn = to N± 7^ 0, to 7^ phase, 
is possibly a weak first order transition. 
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FIG. 2: Left: total gap at CNP for gp = 0.6. Black dots 
are numerical solutions and blue one is the quadratic fit of 
the gap with the field. Right: total gap at CNP for gp = 
O(brown), 0.6(cyan), 0.63(blue), 0.65(red) 

Discussion: To summarize, we show that by develop- 
ing a sizable magnetization in the direction of the field, 
while restricting the Neel order into easy plane one can 
explain the observed scaling of the interaction induced 
gap near CNP in double-gated BLG, at least when the 
field is weak. At moderate field strengths the scaling 
of the interaction induced gap appears to crossover to a 
linear regime, which possibly persists even for stronger 
fields. However, the slope of the linear scaling is off by 
an factor of 1.85 from the experiment. The reason for 
such mismatch is not clear, although it may arise due 
to the existence of additional OPs in high fields. Within 
the zeroth LL, the quantum spin Hall insulator is exactly 
proportional to the magnetization. The electron-electron 
interactions can therefore develop the spin Hall order at 
stronger fields, when the LL quantization is sharp. Yet 
another order, pointed out by Kharitonov '21|, preserving 
the layer inversion symmetry, odd under the exchange of 
Dirac points, and breaks spin rotational symmetry, can 
coexist with the AF order within the zeroth LL. Other- 
wise, it scales linearly with the field and vanishes in its 
absence. Therefore, developing a sizable component of 
this order one can meet the slope at high fields, observed 



experimentally. However, such an order commutes with 
the AF order, unlike the FM order. Therefore existence 
of such order cannot change the curvature of the gap from 
the one with pure AF phase, at weak fields. 
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FIG. 3: Left: Total gap at CNP as a function of the field 
< B{T) < 4.5 for gp ~ 0.65. Red dots are numerically 
found solutions, whereas the solid lines are fitting functions. 
Right: Critical total field (Sf /Bx) for various strength of 
perpendicular field {B±), necessary to destroy the AF order 
at CNP. 

Validity of the proposed scenario can be tested by plac- 
ing the dual-gated BLG in tilted magnetic fields. AF 
order scales with the perpendicular component of the 
field {B±), while the magnetization scales with the to- 
tal magnetic field (Bt). With increasing Bt but a fixed 
B±, the Neel order gradually decreases, whereas the FM 
order increases with Bt and beyond a critical total field 
(JSf ), system can be found in a pure FM phase. With 
gp = 0.65, we here compute B^ for various B±, shown 
in the right column of Fig. [3]. It is evident that for 
stronger B±, B^ gets larger, due to the monotonic in- 
crement of AF order with i3j_. Otherwise, i?f scales 
quadratically with Bj^ , revealing the quadratic variation 
of AF with B±^ , shown in Fig. [1] . With zero field gap 
Ao = 1.13 meV asinRef. [9], Bf « 1.02,2.4,4.8,8.1,15 
T respectively for B^ = 0.68, 1.13, 1.6, 2.03, 2.5 T. In ex- 
periment, one can increase Bt while keeping B±^ fixed by 
tilting the BLG sample, where the tilting angle reads as 
Ot = CQS^'^ {B ^ / Bt) . For a fixed Bj_, B^ is lowered if 
the zero field gap (Aq) gets smaller. For example, with 
Bx_ = 1.6 T, Bf = 4.8, 3.8, 2.95 T for Aq = 1.13, 0.9, 0.7 
meV, respectively, if gp = 0.65. In dual-gated BLG, one 
can reduce the size of the zero field gap by applying a 
weak electric field among the two layers [l^ . Appearance 
of such pure FM phase can be confirmed by measuring 
the longitudinal resistivity R^x [23|- The field induced 
FM phase can host counter propagating edge modes of 
opposite spin polarity, similar to MLG [23], hence R^x 



becomes much smaller than that in presence of a Neel 
order. 

Placing the chemical potential close to the first excited 
state at ±Egap{v = 0), additional incompressible Hall 
states at filling v — ±2 can be formed by developing a 
third component of the Neel order {N-^ = iV||), parallel 
to the applied magnetic field. To the linear order in N^ 
the activation gap for z/ = ±2 reads as 



Egap{v = ±2) = 2 



\ + M 



{nI + (A + MY 



1/2 



iV3 + 0(7V|). 



(14) 
±1 Hall states can also 
[Oj. The scaling of the 



Similarly, in MLG, the v = 
be formed by developing N^ 
V = 1 state in MLG, has been computed for large num- 
ber of fermionic species [2J|. The component of Neel 
order in the direction of the field arises only from the 
zeroth LL, whereas its easy plane component is finite in 
all the LLs. Therefore N^ ~ B, whereas M ^ B and 
{N±{B) — Ao) ~ B^, at least when the field is weak. 
Therefore the gap for v = ±2 Hall state is expected to 
scale quadratically with weak fields, otherwise vanishes 
at zero field. The coefficient of the term proportional 
to B^ in Eq. [14] scales as I/VAq. Therefore curvature 
of the activation gap of i' = ±2 state in weak fields is 
expected to increase upon applying a weak electric field 
between the layers, which decreases the size of the zero 
field gap [1J|. On the other hand, in the presence of 
a strong tilted field, when N± = 0, activation gap for 
1/ = ±2 Hall state scales linearly with the field. A sim- 
ple estimation yields that for weak fields B < 10 T, 
N^ ^ [a/lBY'lJ-B ^ (10^^ — 10~^)/iB per electron, similar 
to that in MLG ll]. Measurement of such a small com- 
ponent of the Neel order &t v = ±2, can test the validity 
of the nature of broken symmetry phases in BLG subject 
to magnetic fields. 
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